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In models of real scalar fields with degenerate double-well potentials, spherically symmetric, 
large amplitude fluctuations away from the vacuum are unstable. Neglecting interactions with 
an external environment, the evolution of such conflgurations may entail the development of an 
oscillon; a localized, non-singular, time-dependent configuration which is extremely long-lived. In 
the present study we investigate numerically how the coupling to a heat bath influences the evolution 
of collapsing bubbles. We show that the existence and lifetime of the oscillon stage is extremely 
sensitive to how strongly the field is coupled to the heat bath. By modeling the coupling through a 
Markovian Langevin equation with viscosity coefficient 7, we find that for 7 ^ 5 x 10~'*m, where m 
is the typical mass scale in the model, oscillons are not observed. 

PACS: ll.lO.Lm, 05.70.Lm, 98.80.Cq 



*NSF Presidential Facuhy Fellow 

t email: gleiserQpeterpan. dartmouth.edu 

*email: rhaas@phys.ufl.edu 



1 



I. INTRODUCTION 

It is well-known that nonlinear field theories allow for 
the existence of static, regular, and localized configura- 
tions [0. In (l+l)-dimensions, one finds the kink so- 
lutions to the Klein-Gordon equation with sine-Gordon 
or 0^ potentials. Unfortunately, for a higher number of 
spatial dimensions, the only static and stable solutions 
involve either two or more fields with some topological 
conservation law, as in the case of the 't Hooft-Polyakov 
monopole or the Nielsen-Olesen vortices or a con- 
served global charge, as in the case of nontopological 
solitons [||. The interest in such static solutions stems 
from a myriad of possible applications, from the modeling 
of particles and cosmological topological defects, to the 
propagation of information in optical fibers and defects 
in liquid crystals and superfiuids. 

Given the justified interest in static solutions, unsta- 
ble, time-dependent solutions have received considerably 
less attention in the literature. However, as recent work 
has emphasized, the presence of nonlinearities can lead 
to the existence of extremely long-lived (i.e., nearly non- 
dissipative) solutions of the Klein-Gordon equation . 
These spherically symmetric solutions, known as oscil- 
lons, are characterized by a rapid oscillation of the field 
at small radial values (the core) of the configuration, very 
much like the (l+l)-dimensional breathers which form 
from kink-antikink collisions In fact, oscillons can 
be thought of as a possible stage during the collapse of 
an unstable bubble; as the bubble collapses radiating its 
initial energy to infinity, it may (or may not) settle into 
an oscillon before completely disappearing. [The detailed 
conditions for a bubble to settle into an oscillon configu- 
ration are given in Ref. 1^.] 

Apart from adding to our knowledge of time-dependent 
coherent nonlinear phenomena in field theories, the fact 
that collapsing bubbles may settle into long-lived config- 
urations may have important consequences for our un- 
derstanding of the dynamics of phase transitions. As an 
example, consider a system in a metastable state which 
is being cooled down at some rate per unit volume, Fcooi- 
For a strong enough first order phase transition, where 
homogeneous nucleation theory is applicable, the system 
will supercool in this metastable state before it decays 
by nucleating a critical radius bubble. The rate per unit 
volume for nucleating a critical bubble is dominated by 
the Boltzmann factor, Feb ^ Gxp[~Fch/T], where Fch is 
the cost in free energy for nucleating the critical bubble, 
and T is the temperature. Critical bubble nucleation is 
suppressed for Fcb/Fcooi < 1- Unless the barrier disap- 
pears below a certain temperature, the free energy cost 
for nucleating a critical bubble typically decreases as the 
system is cooled down, reaching a minimum value before 
it starts increasing again Q. Apart from critical bub- 
bles, smaller free energy, unstable (subcritical) bubbles 



are also nucleated. Since subcritical bubbles may fall into 
an oscillon stage during their collapse, they may exist for 
a much longer time than one would naively estimate. If 
their lifetime, Tosc, is longer than the cooling time-scale 
[tosc > (FcooiV^)"^], it is possible that as the temperature 
decreases, they will overcome the free-energy barrier and 
initiate the decay of the metastable state. In other words, 
oscillons could become critical bubbles. 

In practice, this possibility has not yet been analysed in 
much detail. Estimates in the context of the electroweak 
phase transition show that oscillons are sufficiently sup- 
pressed to be of negligible impact However, based on 
the qualitative arguments above, one expects oscillons to 
be of greater relevance in earlier (in a cosmological con- 
text), and thus faster, phase transitions. In fact, in the 
work of Copeland, Gleiser, and Miiller, it was remarked 
that for a GUT-scale transition oscillons are efficiently 
produced by thermal processes if Fosc/T < 10. This 
condition may be satisfied by, e.g., Coleman- Weinberg 
potentials. Also, apart from their potential relevance to 
the dynamics of both cosmological and laboratory phase 
transitions, it remains to be seen if oscillons exist in the 
non-relativistic limit of the Klein-Gordon equation, as 
time-dependent solutions of the nonlinear Schrodinger 
equation. 

The above discussion neglected the effects of an exter- 
nal environment on the evolution of shrinking bubbles. 
However, in most realistic situations, a thermal back- 
ground will infiuence the dynamics (and coherence) of 
any field configuration. In this work, we will examine 
the evolution of coherent field configurations coupled to 
heat baths in more detail. The evolution of the scalar 
field in a thermal bath will be modelled by a Markovian 
Langevin equation. As we will see, the presence of a 
thermal bath adds additional constraints on the possible 
existence and lifetimes of oscillons; even if the collapsing 
bubble does settle into an oscillon, the bath can affect 
its lifetime quite appreciably. However, for sufficiently 
small couplings between the field and the thermal bath, 
oscillons are still present. 

The rest of this work is organized as follows. In the 
next section, we briefly review the properties of oscillons 
in the absence of a thermal bath. In Section HI, we dis- 
cuss how we study bubble evolution in the presence of 
a heat bath. In Section IV, we present our numerical 
results. Lifetimes of the oscillons for several sets of pa- 
rameters are obtained using a method consistent with an 
ensemble-averaging procedure, together with empirically 
determined equations of best fit for the data. In Section 
V we summarize our results, pointing to possible direc- 
tions of future work. 
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II. COLD OSCILLONS 

In this section we briefly review some of the results 
which established the existence of oscillons as a possible 
stage during bubble collapse for nonlinear field theories. 
This discussion will help us set up the notation and con- 
cepts which will be useful later when we include thermal 
effects. More details can be found in Rcfs. QJ^. 

The action for a real scalar field in (3+l)-dimcnsions 

is 



Figure 1 



(1) 



where we restrict our investigation to a symmetric 
double-well potential of the form 



(2) 
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A solution 0(x, t) to the equation of motion, 



has energy 

E[(t)] ^ f (fx 



(3) 



(4) 



Since we are only interested in spherically-symmetric 
configurations, it proves convenient to introduce dimen- 
sionless variables, p — rm, t — tm, and $ = -^(t)- The 
nonlinear Klein- Gordon equation is. 



92$ 



^2$ 



2(9$ 
P dp 
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(5) 



Oscillons arc found by solving the above equation with 
a bubble-like initial configuration, which so far has been 
taken to be either a "Gaussian" (equ. 6) or a "tanh" 
(equ. 7) bubble, 



*(P,0)= -[($o-*c) tanh(p 



(6) 

i?o) + $o + *c] , (7) 
i?o > 1 . 



i?o is the initial radius of the configuration, $0 is the 
asymptotic vacuum the configuration decays into (we will 
take $0 ~ ^1)5 and $c is the value of the field at the 
configuration's core. [We will take it to be $c = 1, 
the bubble interpolates between the two vacua.] We also 
impose the boundary conditions, 

^ 00, r) = $0, *'(0, r) = 0, 0) = . (8) 



FIG. 1. Bubble energy vs. time for 7 = and T : 
case of (a) R = 2.3, (b) R = 3.4, and (c) R = 2.7. 
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Equation ^ is then numerically solved by a finite dif- 
ference routine second order accurate in time and fourth 
order accurate in space. [We used the same discrete steps 
as in Ref. [^, that is, a spatial step h = 0.1 and a time 
step St = 0.05.] In order to examine the rate at which the 
initial configuration decays into the vacuum, we measure 
the total energy within a spherical shell surrounding it. 
As discussed in Refs. the oscillon configuration is 

identified by the presence of an almost flat plateau in the 
plot of energy vs. time, signalling a regime during which 
almost no energy is radiated away. Typically, lifetimes 
range from lO^m""'^ to lO^m"-'^. In Figure 1, we show the 
energy of a few oscillons as a function of time for several 
Gaussian bubbles of differing initial radii. 



III. HOT OSCILLONS 

We would like to extend the previous analysis for fields 
interacting with thermal baths. The simplest way to do 
this is to couple the scalar field to a thermal bath via a 
generalized Langevin equation, which assumes the bath 
to be Markovian (white noise) and that the field cou- 
ples to the noise additively. The bath is characterized by 
a viscosity coefficient 7 and by a random noise ^ (x, t) , 
which are related by the fluctuation-dissipation theorem. 



(e(x,t)e(x',t')> =27r5(t-t')<5=^(x-x') 



(9) 



7 represents the coupling strength between field and 
thermal bath. For simple models, it can be expressed 
perturbatively in terms of the coupling constants in the 
system. [See, e.g., Ref. [0].] However, we will treat it 
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here as a free parameter. Note that defines the ther- 
mahzation time-scale. More comphcated nonlocal forms 
of the Langevin equation could be used, although we 
think it wise to defer such approaches to later studies. 
[More details can be found in Rcfs. and references 

therein.] 

In terms of the dimcnsionless coordinates defined be- 
fore, the Langevin equation is given by. 



Figure 2a 
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where 7 = 7/TO and ^ = ^/X^/m^ are the dimcnsionless 
viscosity and noise respectively. 

Here, one last simplification comes in. In general, one 
should be solving the fully (3-|-l)-dimensional Langevin 
equation, as the bath carries no symmetries with it. How- 
ever, we will assume that the dominant interactions of the 
bubble with the thermal bath are in the radial direction. 
We must rewrite the fiuctuation-dissipation relation con- 
sistently with spherical symmetry, effectively keeping the 
problem one dimensional. The full (3-|-l)-dimensional 
problem proves to be very CPU intensive and of limited 
interest. [Note that in the presence of a thermal bath, re- 
sults are obtained after an ensemble averaging procedure. 
In other words, each measurement can involve hundreds 
of runs. As it is, the results obtained here already re- 
quired a network of workstations running for a total of 
6000 days of CPU time.] This approximation is not only 
economical, but also makes sense physically. For bubble- 
like configurations, it is reasonable to expect that the 
effect of the bath will be dominant in the radial direc- 
tion. In any case, the results here should be considered 
as an upper bound on the lifetimes; including perturba- 
tions in all three dimensions will not help increase the 
oscillon's lifetime. 

The spherically-symmetric fluctuation-dissipation re- 
lation is. 



{i{p,r)ap',r')) 



1 



2'Kp- 



n9S{T~T')S{p-p') , (11) 



where 9 = X T/m is the dimcnsionless temperature. 

Before proceeding, we must make sure that our expres- 
sion for the fluctuation-dissipation relation is physical. A 
simple way to do this is to consider the evolution of a free 
field in a parabolic potential and measure the energy per 
degree of freedom E/N, where N is the total number of 
spatial lattice points. Since the problem is one dimen- 
sional, the equipartition theorem states that in equilib- 
rium, E/N = r/2. We have confirmed that this relation 
is satisfied for a wide range of parameters, as shown in 
Figure 2. 
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Computed energy per degree of freedom, tc 



7 ^ 



10" 



and r = 10.0, (b) 7 = 0.1 and T = 0.1, 
(c) 7 = 1.0 and T = 0.25. The horizontal lines in each figure 
represent the expected values from the equipartition theorem. 
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IV. NUMERICAL RESULTS 

Now that we have an evolution equation in the pres- 
ence of a thermal bath consistent with the assumption of 
spherical symmetry, we can start investigating how the 
bath influences the formation and longevity of oscillons. 

As in the "cold oscillon" case, we measure the energy 
within a shell surrounding the initial configuration as a 
function of time. The radius of the shell, i?shoU, should 
be taken to be sufficiently larger than the radius of the 
initial configuration, Rq. Typically, we chose i?shoii = 15. 
The presence of temperature, however, immediately in- 
troduces two difficulties. First, the bath contributes to 
the energy within the shell, obscuring the measurement 
of the bubble's energy. Second, a physical result must 
be obtained by an ensemble averaging procedure, where 
different conditions are set by different seeds for the ran- 
dom number generator responsible for the noise term in 
the equation. However, as soon as we change the seed, 
the oscillon lifetime changes, as we show in Figure 3 for a 
selection of 30 runs out of 150. We must find a consistent 
procedure for measuring the ensemble-averaged oscillon 
lifetime. Let us deal with each of these problems in turn. 



Figure 3 
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FIG. 3. Bubble energy vs. time for R = 3.2, 7 = 10"^, and 
T = 0.1 with different seed values for each curve. The figure 
represents only 30 runs from a total of 105. 

The first problem can be taken care of quite easily; 
we perform several runs without a bubble [i.e., with a 
constant initial configuration, <i>(/9, 0) = — 1], and mea- 
sure the ensemble-averaged energy of the bath within the 
shell, (-Ebath)- The angular brackets denote ensemble- 
averaged quantities. We then define the normalized en- 
ergy of the bubble by 

Ebuh = i^total — (Ehath) , (12) 



where, Stotai is the total energy within the shell. Figure 
4 shows a typical case for E'bub, -Etotai, and (i^bath)- As a 
test of this procedure, we confirmed that (-Ebub) — > 
after bubbles have decayed. Another test is to mea- 
sure the plateau energy for small 7 and T. The result 
should be (and is) consistent with the cold oscillon case, 
E^piatoau ~ 43to/ a. 



Figure 4 
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FIG. 4. Bubble energy, total energy, and bath energy as 
functions of time for Ro = 2.7, 7 — 10~^, and T — 5.0. 

In order to handle the second problem, we first define 
how we measure the oscillon's lifetime. As opposed to 
the case of cold oscillons, where we could just estimate 
the lifetime by checking when the oscillon finally decays 
(see Figure 3), here we must automate the procedure, as 
we will have to handle thousands of runs. First, we must 
establish if a given initial bubble settles into an oscillon 
stage. As we remarked earlier, the most important sig- 
nature of the oscillon is the near constancy of its total 
energy. In Refs. it was measured that during the 

oscillon stage the plateau energy varies by at most 20%. 
After an exhaustive search, we found that we can identify 
an oscillon stage by setting Id-Epiateau/c?*! < 0.02m'^/A. 
Steeper changes of the plateau energy were incompatible 
with the nearly non-dissipative nature of oscillons. Now 
that we know how to hunt for oscillons, we must de- 
cide how we will measure their lifetimes. Here, the fact 
that these configurations decay rather quickly allows us 
to simply set the oscillon lifetime to the time when its 
energy drops to O.liJpiatoau- 

To summarize so far, we first choose values for the three 
variables in the model: the initial radius Rq , the viscosity 
7, and the temperature T. Each numerical "experiment" 
consists of evolving the Langevin equation with a choice 
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of variables and a given value for the random number 
seed. The code then searches for an oscillon by identify- 
ing a fairly flat plateau, and then measures its lifetime by 
recording the time when its energy decays to O.liSpiatoau- 
This procedure is repeated a large number of times with 
different seeds for the random number generator. This 
way, we obtain a distribution of lifetimes for a given set 
of variables, as shown in Figure 3. In order to determine 
the ensemble-averaged lifetime, we produce a histogram 
where we bin different decay times vs. the number of 
occurrences. In Figure 5, we show a few examples of 
histograms for varying radii and fixed viscosity and tem- 
perature. In Figure 6, we show a few examples of varying 
viscosity, with fixed radius and temperature. 
Figure 5a 



Figure 6a 




Figure 5b 





Figure 6b 




FIG. 5. Histograms for (a) R = 2.5 and (b) R = 3.0, for 
7 = 10~^ and T = 0.1. At indicates the bin size used in the 
specific cases. 



FIG. 6. Histograms for (a) 7 = lO"'^ and (b) 7 = 3 x 10"", 
for Rf) = 2.7 and T = 0.1. At indicates the bin size used in 
the specific cases. 



We can now extract the ensemble-averaged lifetime 
by fitting a Gaussian to the various histograms, using 
a least-squares approximation. Since the thermal noise 
producing the variation in the lifetimes is generated from 
normal deviates, it is reasonable to assume that the life- 
times will follow the central limit theorem, and produce 
Gaussian distributed data as well, for a sufficiently large 
number of runs. The mean gives the desired lifetime and 
the variance gives the spread in the lifetime measure- 
ments. In order to obtain the best possible fit, we used 
the Levenberg-Marquardt method 

The results for the lifetime as a function of radius and 
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different values of viscosity and temperature are shown 
in Figure 7. It is clear that the coupling to the bath 
strongly suppresses the duration of the oscillon stage, 
when it is at all present; for values of viscosity larger 
than 7 > 5 X 10~^m, no oscillons could be found. Bub- 
bles quickly collapse, radiating their coherent energy into 
the incoherent thermal bath. Since is the typical 
thermalization time-scale, we obtain a rough criterion for 
the existence of oscillons in the presence of a Markovian 
thermal bath, 



7 > ^0 



(13) 



where tq is the lifetime of the cold oscillon of the same 
radius. Since typically 10'^ < tq < 10^, this result is con- 
sistent with the fact that the bath acts to destroy the 
coherence of the field configuration. 

Figure 7 




FIG. 7. Oscillon lifetime vs. radius for constant viscosity 
and temperature. 

From Figure 7 we note that although the bath affects 
the longevity of the oscillons, it does not appreciably af- 
fect the range of radii for the initial bubble settling into 
the oscillon stage; as in Refs. 0,^, we still obtain roughly 
2.3 ^ Rq 4.5 for the allowed range. Also, the longest 
lived oscillons are the most affected, consistent with the 
above condition, Eq. 

In Figure 8 we display the results for lifetime vs. vis- 
cosity, for Rq = 2.7m~^ and T = 0.1m/ X. A least- 
squares fit produces the empirical relation. 



t(7) = A 



1 - 



B 



70 



-c 



(14) 



where A = 2.0 x 10^, B = 4.0 x 10"^ C = 0.43, and, 
7o = 5.8 X 10~^. The error bars are given by the variance 



of the Gaussian fit the histograms. 70 sets the scale for 
the existence of oscillons. 

Figure 8 
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FIG. 8. Oscillon lifetime vs. viscosity for constant ra- 
dius and temperature. The equation providing best fit is 
A (^)~^ where yl = 2.0 x 10^ B = 4.0 x 10"^ C = 0.43, 
and 70 = 5.8 x 10"^ 



Figure 9 
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FIG. 9. Oscillon lifetime vs. temperature for constant 
radius and viscosity. The equation providing best fit is 
^exp {-jr) + B where A = 2.1 x 10^, B = 1.7 x 10^, and 
To = 2.1. " 



Finally, in Figure 9 we display the results of lifetime vs. 
temperature, for Rq = 2.7m~^ and 7 = 10~^. It is clear 
that the lifetime is less sensitive to temperature then to 
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viscosity. Again, an empirical fit can be obtained, 

r(r) =Aexp(^-|r) +S , (15) 

where A = 2.1 x 10^, B = 1.7 x 10^, and Tq = 2.1. For 
higher temperatures, it becomes very difficult to extract 
meaningful results, as the signal to noise ratio sharply 
decreases. 

V. CONCLUSIONS 

We investigated the evolution of collapsing bubble-like 
configurations in the presence of a thermal bath. The 
bubbles were assumed to have Gaussian profile, and the 
dynamics were modelled by a generalized Langevin equa- 
tion with a Markovian thermal bath. In the absence of a 
thermal bath, it is possible for these configurations to set- 
tle into what is known as an oscillon stage, a nearly non- 
dissipative configuration which is thus extremely long- 
lived. 

We showed that the presence of a thermal bath af- 
fects the existence and longevity of the oscillon stage 
rather strongly. For values of the viscosity coefficient, 
7 > 5 X 10~^TO, no oscillon stage develops, independent 
of temperature. This is naively analogous to an over- 
damped harmonic oscillator. For smaller values of 7, 
an oscillon stage develops, but its lifetime is suppressed. 
This result can be related to the fact that the time-scale 
for thermalization is ~ 7~^, while in the absence of a 
heat bath the typical lifetime of the oscillon stage is 
To ~ 10'^ — 10'^m~^. This is analogous to an under- 
damped harmonic oscillator, where the decay time-scale 
is 7^^. When the decay time-scale is of order the life- 
time of the oscillon, the bath suppresses the duration of 
the oscillon stage. For even smaller values of 7, the os- 
cillon approaches similar lifetimes as in the cold oscillon 
case. Thus, the relevance of oscillons in different con- 
texts in which a thermal bath is present will depend on 
the strength of the coupling between the field and the 
thermal bath, here modelled by 7. As a rough rule, os- 
cillons are present whenever 7"^ > tq. 

Our discussion focused on Markovian thermal baths. 
However, given the nature of interacting field theories, 
it is possible that the bath will have more complicated 
properties such as spatio-temporal correlations which are 
nonlocal, and/or couplings which are multiplicative as 
oppose to additive. (For example, terms of the form ^0 
in the equation of motion |^.) It remains to be seen 
if these nonlocal effects will act as an additional source 



of coherence for fields. That being the case, coupling 
to more general thermal baths may have very different 
consequences for the existence and longevity of oscillons. 
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